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Abstract. Derjugin’s Mountain Newt (Neurergus derjugini derjugini) inhabits the Zagros Mountains in northeastern Iraq 
and western Iran. The present study was conducted in Rikhalan village (Mirisour stream) in Marivan County, western 
Kurdistan Province, to determine habitat requirements and population size of Derjugin’s Mountain Newt. From April 
through May 2020, the capture-mark-release-recapture technique was used to estimate population size. Habitat param- 
eters including pH, dissolved oxygen and electrical conductivity of water, water and air temperature, the size of bedrocks, 
water flow speed, and canopy and aquatic vegetation cover were measured. A total of 3175 specimens of Derjugin’s Moun- 
tain Newt (1853 males, 1322 females) were captured and recaptured. Population size was estimated using the software CAP- 
TURE and the Jolley-Seber method. According to the results obtained by means of the MARK and Ecological Methodol- 
ogy software, the population size of Derjugin’s Mountain Newt in the study area was between 2401 and 2554/1664 m7”, 1165 
and 2877/1664 m°, respectively. Habitat preferences of Derjugin’s Mountain Newt in the Mirisour Stream are clear shallow 
water with a slow flow speed, pH 8.23-8.36, electrical conductivity 292-331 uS/cm, dissolved oxygen 7.26-8.42 mg/l, water 
temperature 12-17.5°C, and vegetation to include oak (Quercus brantii, Q. infectoria), plane (Platanus orientalis), walnut 
(Juglans regia), and willow (Salix spp.) while the aquatic vegetation consists of broadleaf (Mentha sp.) and orchard grass 
(Dactylis glomerata). Presence and abundance of this species were considerably higher at greater altitudes. The key threat 
to the species is habitat loss, which is caused by water being siphoned off for agricultural purposes and blocking springs 
as a result. 
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habitat variables. 


Introduction 


Amphibians are an integral part of many natural ecosys- 
tems. They serve as main prey for predators such as birds 
and reptiles across different ecosystems (WOOLBRIGHT 
1991, DUELLMAN & TRUEB 1994). Decreasing amphibian 
populations could have a major impact on natural eco- 
systems (GARDNER 2001). Despite their ecological values, 
various anthropogenic factors currently affect their popu- 
lations either directly or indirectly, and their populations 
are declining in many parts of the world. Amphibians, 
as much as other animal species, suffer from habitat de- 
struction and fragmentation globally (ASHRAFZADEH et 
al. 2019). Derjugin’s Mountain Newt, Neurergus derjugini 
(NESTEROV, 1916) inhabits the western margin of the Za- 
gros Mountains in Iran and northeastern Iraq in a narrow 
distribution range. This species has two subspecies: Der- 
jugin’s Mountain Newt (Neurergus d. derjugini) and Small- 
spotted Mountain Newt (Neurergus d. microspilotus) (Sa- 


FAEI-MAHROO & GHAFFARI 2020). THE IUCN Red List of 
Threatened Species classifies this species as ‘critically en- 
dangered’ (SAFAEI-MAHROO et al. 2015). 

A key parameter for evaluating the status of species con- 
servation and the extinction risk level is population size 
(SEGEV et al. 2010). Population monitoring is one of the 
principles of wildlife management. Undoubtedly, the lack 
of information about population size, abundance and den- 
sity of species impedes informed conservation manage- 
ment (KREBS 1989). Also, habitat variables need to be con- 
sidered in management plans (UNGLAUB 2015). Population 
size (i.e., Minimum Viable Population) and habitat require- 
ments are vital to safeguard the survival of endangered spe- 
cies. In the present study, the population size of Derjugin’s 
Mountain Newt (Neurergus d. derjugini) was estimated 
based on the capture-mark-release-recapture method, 
and some variables in the aquatic and adjacent terrestrial 
habitats in Rikhalan village (Mirisour Stream) in Marivan 
County, western Kurdistan Province, were investigated. 
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Material and methods 
Study site 


This study was carried out in a stream called Mirisour, 
Marivan County, Kurdistan Province, in western Iran 
(35°26.438” N, 46°9.2’68” E) (Fig. 1). The climatic con- 
ditions of the study area are categorized as cold and hu- 
mid with an average temperature of 12.8°C and the aver- 
age annual precipitation amounts to 991.2 mm (HANAFI & 
HATAMI 2013). The study site is located at an altitude of 
1320-1555 m above sea level where the vegetation includes 
oak (Quercus brantii, Quercus infectoria), walnut (Juglans 
regia), and willows (Salix spp.). 


Population estimation 


Fieldwork was conducted from 22 April through 17 May 
2020. To estimate the population size, capture-recapture 
sampling was performed once every five days by four per- 
sons in six surveys in which adult specimens were captured 
only. Population estimation was performed by means of the 
Jolly-Seber method in Ecological Methodology software 
(version 7.2). The study area was divided into three sections 
with lengths of 170 (the first station), 270 (the second sta- 
tion) and 600 metres (the third station), respectively. Sam- 
pling was always conducted by walking along the stream 
for 3 hours (from 11 through 14 h) (SAEAEI-MAHROO et al. 
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Figure 1. Location of the study area. 
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2020). Specimens were captured by hand and were marked 
by toe clipping before being released in their original cap- 
ture locations. During the first visit, the third toe of the left 
foot, in the second visit, the third toe of the right foot, in 
the third visit, the second toe of the left hand, in the fourth 
visit, the second toe of the right hand, and in the fifth visit, 
the first toe of the left foot were cut off. The numbers of new 
captures and recaptures were recorded on each occasion 
(Table 1). In Table 1, only recaptures with one mark were 
recorded in the recapture column and samples with two 
or more marks were considered as total recaptures. There- 
fore, to estimate the population size of Derjugin’s Mountain 
Newt, those estimates were considered that had the high- 
est ratios of markers and the lowest standard errors. Ac- 
cordingly, the estimated population values were tabulated 
(Table 2) and then each station was considered a replica- 
tion, for which averaging was performed. To estimate the 
total population size of adult newts, a closed population 
was designated. We then applied the software CAPTURE in 
MARK (version 9.0) (Oris et al. 1978) and the Jolley-Seber 
method in Ecological Methodology software. To choose the 
best model for estimating population size, the Goodness of 
fit method in MARK software was used (O’DONNELL et al. 
2008). As each of these population estimation models has 
an error if the value of p_ hat (probability of capture) is 
equal to or greater than 0.1 and the value of t (number of 
trapping occasions) is equal to or above 5 (p 2 0.1, t 25), the 
bias due to population size estimation by the models can be 
neglected (Oris et al. 1978). According to Table 3, the value 
of p_ hat is equal to 0.2, and t is equal to 6. Selecting and 
testing the appropriate model for capture and recapture 
methods in closed populations is still an unresolved issue 
(CoocH & WHITE 2014), though. If only one of the varia- 
bles, such as time, behaviour or heterogeneity, is considered 
to be effective in capture probability, then Goodness of fit 
can be considered a suitable model for testing the Mo, Mt, 
and Mb models. As suggested by Table 4, model Mo vs. Mt, 
means that time is the factor of change in the probability of 
catching and model Mo vs. Mb means the behaviour of the 
factor of change in the probability of capture. 


Habitat variables 


As a part of the study, the stream over a length of 1109 m 
was selected as a macrohabitat to measure the variables of 
aquatic habitat. Then, this section was divided into three 
consecutive stations as microhabitats to determine the var- 
iations between habitat characteristics and the relationship 
between species density and habitat characteristics. In our 
field surveys, habitat variables including water temperature 
and air temperature, measured by glass thermometer, elec- 
trical conductivity, acidity and dissolved oxygen, through 
the multi-parameter device (model: HACH, HQ4od), alti- 
tude, stream area and length and latitude, by the GPS de- 
vice (model: GARMIN GPS MAP 788), slope percentage 
(Compass software) and water flow velocity were meas- 
ured and recorded. To determine the water flow veloc- 
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Table 1. Counts of new captures and recaptures in the first to the sixth samplings. 


Station Date Total captures New captures  Recapture 1 
1 22/4/2020 72 - - 
2 22/4/2020 168 - - 
3 22/4/2020 382 - - 
1 27/4/2020 77 46 31 
2 27/4/2020 135 77 58 
3 27/4/2020 227 188 36 
1 2/5/2020 73 34 16 
2 2/5/2020 115 50 30 
3 2/5/2020 343 225 70 
1 7/5/2020 63 24 10 
2 7/5/2020 136 50 30 
3 7/5/2020 294 183 37 
1 12/5/2020 44 12 2 
2 12/5/2020 118 24 19 
3 12/5/2020 339 161 31 
1 17/5/2020 27 7 2 
2 17/5/2020 105 18 11 
3 17/5/2020 443 169 32 


Recapture2 Recapture3 Recapture4 Recapture 5 


11 = - z 
19 : - - 
38 = = = 
4 10 - = 
20 8 7 7 
34 19 - = 
= 6 6 = 
13 11 16 = 
26 40 24 = 
2 - 2 2 
6 7 6 6 
29 33 35 39 


Table 2. Population size at each station as estimated with Ecologi- 
cal Methodology software. 


Proportion 


; F A 5 
Station Estimate+SE marked Lower 95% Upper 95% 
389+91.3 0.41 269 715 
1129+455.7 0.5 628 2822 
3393+575.4 0.35 2600 5095 


Table 3. Population size estimate for Derjugin’s Mountain Newt 
in the study area with a total area of 1664 m°. t = number of 
trapping occasions, Mt+1 = number of animals captured, P-hat 
= estimated probability of capture, CI = an approximate 95% 
confidence interval. 


Model Estimator t Mt+1 P-hat Estimate+SE 95% CI 

MO - 6 1892 0.2136 24774+39.38 2406-2560 
Mt Chao 6 1892 0.2066 25364+54.86 2437-2652 
Mt Darroch 6 1892 0.2133 2472+39.10 2401-2554 
Mb - 6 1892 0.2192 24454+66.68 2330-2592 


Table 4. The model selected by the Goodness of fit method. 


M0 vs. Mh MOvs.Mt MO vs. Mb 
Chi_ Square 5.47 67.4 0.3 
Degrees of freedom 3 5 1 
Sig 0.14 0.000 0.56 


Table 5. Criteria for model selection by the Goodness of fit 
method. 


Mt Mb 
1 0.13 


Model criteria 


ity, two points between a certain distance were selected. 
After measuring the distance between these two points 
(in m), a piece of wood was dropped into the water at the 
upstream point and the time required for it to reach the 
downstream point was taken with a stopwatch. Finally, the 
measured distance was divided by the time recorded and 
water flow velocity was calculated based on m/s (ALLAN 
2007). It should be noted that the velocity of water flow 
was measured at several points and extrapolated to apply to 
the whole stream. Water depth was also measured at vari- 
ous points by a metal measuring tape. Also, in each station, 
the type and percentage of vegetation inside and outside 
the water, the percentage of the canopy cover was recorded, 
and vertebrate (SAFAEI-MAHROO & GHAFFARI 2020) and 
invertebrate (CLIFFORD 1991) species coexisting with Der- 
jugin’s Mountain Newt were noted. The bedrocks of the wa- 
ter course were categorized based on the method suggest- 
ed by SADEGHIFAR & AZARMSA (2016), i.e., bedrocks were 
classified into four different categories: pebble (4-64 mm 
in grain size), cobble (64-256 mm), boulder (> 256 mm), 
and sand-clay (particles < 2 mm). To evaluate the differ- 
ences between physical and chemical factors of the water 
and the difference in altitude between the three stations, 
one-way ANOVA, and to investigate the relationship be- 
tween habitat variables and species density at each station, 
a Pearson correlation test (Pearson’s correlation coeffi- 
cient) were performed using SPSS software (Version 16.0). 
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Table 6. Results of the amount of environmental, physical and chemical factors at each station by one-way ANOVA. 


Station Water Air DO pH EC Altitude 
temperature temperature 
13.9+1.7 20.3+1.8 7.89+0.28 8.31+0.8 296.9410 1335411 
2 13.3+1.5 19.9+1.7 8.44+0 8.3740.8 309.6412 1385412 
3 13.8+1.7 19.641.5 7.5640.18 8.2040.1 336.4414 1554418 
Results 0.98) with density, but this relationship was not significant 


Population estimation 


During our field surveys, a total of 3175 adult Derjugin’s 
Mountain Newts were captured (Table 1). Based on the ob- 
tained results, with a 95% confidence interval, its popula- 
tion was estimated to be between 1165 and 2877 individuals 
in the study area of 1664 m°. Furthermore, the population 
size was estimated with the MARK software (Mt, Mb and 
Mo estimators in CAPTURE) (Table 3) at 2500/1664 m° 
in the study area. Based on the obtained results, Mt (Dar- 
roch) model turned out to be the best model for estimating 
population size (Table 4). According to the Mt (Darroch) 
model, the population size of Derjugin’s Mountain Newt in 
the study area with its area of 1664 m° was estimated with 
a 95% confidence interval to be 2401-2554 individuals, and 
the density was 0.16 to 0.17 individuals/m7. The total area 
of the stream was 14,486 m°. Given the fact that the spe- 
cies was observed only at altitudes above 1370 m, this esti- 
mate could be extrapolated to the entire water area. Thus, 
according to the Ecological Methodology, the population 
size of Derjugin’s Mountain Newt in the entire water body 
with its area of 14486 m? was 1165 to 2877 and the density 
was 0.08 to 0.2 individuals/m?. According to these results, 
Mt is the only model that is significant compared to Mo 
(P < 0.05), and it may therefore be concluded that time is 
the only source of change in the probability of capture and 
Mt can be used as a suitable model for population estima- 
tion. The results of the Goodness of fit model (Table 5) also 
indicate that the Mt (Darroch) model is the best model for 
estimating population size, because it has the highest value 
compared to the other two models. 


Habitat Requirements 


The results of a one-way ANOVA indicated that among the 
water physical and chemical factors, acidity and electrical 
conductivity differed significantly between the three sta- 
tions. There is also a notable difference between stations in 
terms of altitude (Table 6). 

According to the results of the Pearson correlation test, 
the correlation coefficient between acidity and dissolved 
oxygen and the density of specimens was negative (pH = 
-0.89, DO = -0.98). The relationship between pH and the 
density was not significant at a level of 0.05 (P > 0.05, Sig = 
0.3). The relationship between DO and density was not sig- 
nificant at a level of 0.05 (P > 0.05, Sig = 0.1). Electrical 
conductivity had a positive correlation coefficient (EC = 
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at a 0.05 level (P > 0.05, Sig = 0.1). The correlation coeffi- 
cient between altitude and density was positive (= 1), that 
is, density increases with increasing altitude, and the rela- 
tionship between these two variables was significant at a 
level of 0.05 (P < 0.05, Sig = 0.009). Therefore, the values of 
acidity, electrical conductivity, dissolved oxygen, and water 
temperature at the three studied stations are slightly dif- 
ferent from each other yet have no effect on the density of 
samples in the selectedsectors. The altitude variable, which 
was increasing from Station 1 to Station 3, proved to be an 
important factor determining the density of the Derjugin’s 
Mountain Newt in its habitat in the MirisourStream. 
Water flow velocity was measured to be between 0.2 and 
0.5 m/s, and water depth veried between 2 and 67 cm. The 
stream’s bedrock was composed of 40% pebble (4-64 mm), 
25% cobble (64-256 mm), 15% boulder (> 256 mm), and 
20% sand-clay (particles < 2 mm). It was found that the 
species occurred in higher density in areas where the pro- 
portional contents of pebble and sand-clay were larger. The 
dominant vegetation in the study area was oak (Quercus 
brantii, Q. infectoria), plane (Platanus orientalis), walnut 
(Juglans regia), and willows (Salix spp.), and aquatic plants 
were broadleaf (Mentha sp.) and orchard grass (Dactylis 
glomerata). In some parts of the water course, bryophytes 
were also observed (Fig. 2). Aquatic vegetation was sparse 
(30%) and vegetation debris covered a significant percent- 
age of the streambed. Due to the dark brown colour pattern 
of Derjugin’s Mountain Newt larvae, brown accumulations 
of leaf litter are a good place to find refuge. This species 
spawns on submerged branches (Fig. 3). The species coex- 
isting in the same habitat included Bedriaga’s marsh frog 
(Pelophylax bedriagae), freshwater crab (Potamon sp.), and 
insect larvae of the families Gyrinidae, Geriidae Hepta- 
geniidae, Stratiomyidae, Tipulidae and Zygoptera. 


Discussion 


This study and previous research (MIRANI 2013) have found 
that the best time to study the population size of adult 
Derjugin’s Mountain Newt is in the reproduction season 
from early April through late May, a period during which 
both males and females are present in the water. Accord- 
ing to studies conducted by Fr1z1 & EzzaTI (2019), who 
focused on the habitat status and population of the species 
in the Mirisour Stream, the population size estimated by 
the Lincoln-Petersen method was about 1000 in an area 
of 1120 m°. Our results indicate that the topographic char- 
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acteristics of the terrestrial habitat of Derjugin’s Mountain restrial habitats adjacent to the streams ranged from tree- 
Newt around the stream are a steeply sloped forest habitat less floating meadows to dense oak forests. In the study 
with a sparse covering of oak trees in grassland. This spe- area, the values of pH ranged from 8.36 to 8.23, electrical 
cies was present at altitudes above 1,000 m, and its density 
and presence increased with increasing altitude. Based on 
the studies of AFROOSHEH et al. (2016), which examined 
the distribution of the Yellow-spotted Newt (Neurergus d. 
microspilotus) in western Iran, all localities of this taxon 
ranged between 630 and 2057 m above sea level. The ter- 


Figure 3. Camouflaged larvae and spawn on submerged branches, 
17 May 2020. 


Figure 2. Aquatic and terrestrial vegetation: the upper image 
shows aquatic plants, 22 April 2020, and the lower image shows Figure 4. A water pipe draining the stream source for agricultural 
trees around the stream, 17 May 2020. irrigation. 
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conductivity was somewhere between 292 and 331 uS/cm, 
dissolved oxygen registered at 7.26-8.42 mg/l, and water 
temperature varied between 12.5 and 17.5°C. According to 
SHARIFI & ASSADIAN (2004), who examined the distribu- 
tion and conservation status of the Yellow-spotted Newt in 
a stream in Qory Qala in Kermanshah Province, electri- 
cal conductivity was between 323 and 356 uS/cm, pH was 
at 6.65-8.15, and water temperature was reported to be 
11°C. Lowe et al. (2012) examined the physical and chemi- 
cal factors of the aquatic habitat of Gyrinophilus porphy- 
riticus in North America, where the pH was reported to 
be 6-5, electrical conductivity ranged between 12-15 uS/ 
cm, and temperature was between 13 and 17°C. The bed of 
the studied aquatic habitat was mainly composed of peb- 
ble and cobble with 1.5-250 mm in grain size. SHARIFI & 
ASSADIAN (2005) examined the bed of the Kavat River in 
Kermanshah Province, which had the highest number of 
Yellow-spotted Newts. They recorded the size of its cobble 
as 2-256 mm. Our results indicate that Derjugin’s Moun- 
tain Newt occurs at its highest density at higher altitudes 
in the study area. Habitat preferences of Derjugin’s Moun- 
tain Newt in the Mirisour Stream are clear shallow water 
with a slow flow and good oxygenation, pH 8.23-8.36, elec- 
trical conductivity 292-331 uS/cm, dissolved oxygen 7.26- 
8.42 mg/l, water temperature 12-17.5°C, sparse aquatic veg- 
etation, bed material with a high content of sand-clay and 
low content of cobble, and a dense canopy cover. The spe- 
cies is more abundant on the margins and floors of shallow 
areas, where the water flow is slow. Based on SAFAEI- MAH- 
ROO & GHAFFARI (2020) and our observations, the main 
threat to the species is habitat loss imposed by water be- 
ing siphoned off for agricultural irrigation in a manner that 
has effectively blocked the springs (Fig. 4). 
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